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Neutrinos from core collapse supernovae can be emitted from a rapidly accreting disk surrounding 
a black hole, instead of the canonical proto-neutron star. For Galactic events, detector count rates 
are considerable and in fact can be in the thousands for Super-Kamiokande. The rate of occurrence 
of these accreting disks in the Galaxy is predicted to be on the order of ~ 10 _5 yr _1 , yet there is 
little observational evidence to provide an upper limit on their formation rate. It would therefore 
be useful to discriminate between neutrinos which have been produced in a proto-neutron star and 
those which have been produced accretion disks. In order to distinguish between the two scenarios, 
either the time profile of the neutrino luminosity or the relative fluxes of different neutrino flavors 
may be considered. There are some signals that would clearly point to one scenario or the other. 

PACS numbers: 26.50.+x,13.15.+g,97.60.Bw 



Recent advances in the theory of gamma ray bursts 
(GRBs) have pointed to rapidly accreting disks surround- 
ing black holes (AD-BH) because this configuration is 
a natural one for producing massive energy generation. 
Simulations of neutron star mergers and core collapse 
supernovae have indeed found that under certain con- 
ditions, these disks can form 

13 Q Spectra from 
X-ray afterglows of long duration bursts have indicated 
that the bursts are associated with Type Ib/Ic super- 
novae, suggesting massive stellar core collapse as the ori- 
gin. The rate of gamma ray bursts in the galactic region 
is estimated to be approximately 10~ 5 yr _1 with approx- 
imately 2/3 of these being long duration events [4]. This 
is roughly 0.1% of the core collapse supernova rate. Disks 
which produce enough energy for a burst must have fairly 
high accretion rates in the range of 0.1 Mq/s to 1.0 
M Q / s. Observations of short bursts suggest that they 
are associated with compact object mergers, neutron star 
- neutron star, or black hole- neutron star [5|. 

Hypernovae are very energetic core collapse supernova. 
The rate of Type Ic supernovae in an average galaxy 
is about 10 _3 yr _1 , whereas the inferred hypernova rate 
is about 10~ 5 yr -1 Q| . A limited data set of Type Ic 
hypernovae and supernovae shows that supernovae which 
come from progenitors above ~ 20Mq show evidence for 
two branches with a few stars at high energy ~ 10 52 
ergs, and one at low energy ~ 10 50 ergs 6]. One can 
speculate that the more luminous events are produced 
by an accretion disk and the less luminous event is not. 

It would be useful to translate these limits on hyper- 
novae and gamma ray bursts into limits on the AD-BH 
formation rate. Not every AD-BH may make a hyper- 
nova or GRB. An observationally based upper limit can 
be made from the ratio of black holes to neutron stars. 
In the case of low-mass x-ray transients the observations 
will support turning every star above ~ 20M Q into a 
black hole Q, producing a rate of 10% or even greater. 



In this paper we point out that the next Galactic su- 
pernova may produce neutrinos in an accretion disk, not 
in a proto-neutron star. Even if the supernova light curve 
is observed, it would be difficult to use this to determine 
whether a PNS or an AD-BH is formed at the core. In 
this paper we show that the emitted neutrinos may be 
used to immediately and directly make this distinction. 

Detectors that will observe neutrinos from the next 
galactic supernova are currently Super-Kamiokande 
KamLAND_P , MiniBooNE [l(J, SNO [ll| and the future 
Borexino |l2j . A proto-neutron star event will produce 
1000s of events in SuperK, and 100s in KamLAND and 
MiniBooNE. At least one count from an AD-BH super- 
nova may occur in SuperKamiokande even if the event is 
3 Mpc away [13. 

To estimate the neutrino signal, we first examine the 
disks from which they are emitted. Neutrinos in accre- 
tion disks are produced primarily in inverse beta decay 
events. For sufficiently rapidly accreting disks the cen- 
ter of the disk has high enough temperature that the 
neutrinos become trapped. Fig. ^ shows the surfaces 
of last scattering for the neutrinos in the PNS and the 
AD. Fig^i is based on the calculations of |ljj for an AD 
of M = I.OMq/ s, and black hole spin parameter a= 0, 
while Fig. ^ was estimated from the protoneutron star 
temperature and density profile of Ref . |l5j . For a disk 
which accretes at a rate 0.1 Mq/s, the neutrinos are only 
barely trapped. 

If one wishes to detect neutrinos from an AD-BH su- 
pernova, the total energy emitted in the form of neutri- 
nos is clearly important. For the PNS, analyses are often 
done by taking 3 x 10 53 erg of gravitational binding energy 
liberated in the collapse and dividing it equally between 
six species, three flavors of neutrinos and three of an- 
tineutrino. In the case of the AD-BH, gravitational bind- 
ing energy is released as material spirals in from the outer 
edges of the disk but only some of that is released in the 
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FIG. 1: Shows the neutrino surfaces for a proto-neutron star and a 1 Mq/s accretion disk. 



form of neutrinos. The gravitational binding energy re- 
leased, given a 3 M Q black hole, is 9 (M/M ) x 10 53 ergs, 
where M is the mass of material processed by the disk. 
When the neutrinos are not trapped the rate of energy 
loss is determined by the rate at which neutrinos can be 
produced through inverse beta processes. When neutri- 
nos are trapped it is determined by the time it takes the 
neutrinos to scatter out of the trapped region. We find 
that roughly 20% of this gravitational binding energy is 
released in the form of neutrinos for disks of 1.0 Mq/ s 
and 5% in the case of 0.1 Mq/ s. 

The timescale for neutrino emission through the proto- 
neutron star is about 10 s, with an exponentially decaying 
luminosity so that most of the neutrinos are emitted in 
the first three seconds. If an accretion disk is steady 
state, the neutrino luminosity will not fall off but will be 
approximately constant. For example, for a 0.01 Mq / s 
disk which processes 1 Mq, the neutrino emission will 
be spread over 100 s, although deviations from steady 
state will likely occur. The observed time profile of the 
neutrinos may be used as a constraint. 

The energy spectra of the emitted neutrinos will de- 
termine the flux averaged cross sections in supernova 
neutrino detectors. Over the trapping surface in the 
M = 1Mq/s AD, the energies are in the range of 
(E Uc ) ~ 13MeV to (K.) ~ 14MeV and {E„ a ) ~ 16MeV 
to (E^) ~ 18MeV 14] . There is relatively little emis- 
sion due to )jl and r type neutrinos in such disks although 
for faster accretion rates such as 10 Mq/s, fi and r neu- 
trinos are not only produced but also become trapped. 
In more slowly rotating, optically thin disks, the aver- 
age energy of the neutrinos comes close to tracking the 
average thermal energy of the electrons and positrons. 

As can be seen from Fig. [21 the spectra of the v e and v e 
are similar in the PNS and AD-BH cases. The muon and 
tau types however, are quite different, suggesting that a 
detection of the neutrinos by flavor will give a clear sig- 
nature for one scenario or the other. However, neutrino 
flavor transformation must be taken into account when 
translating the flavors of emitted spectra to the flavors 



that will be measured at the detector. 

In both the PNS and AD-BH, the neutrinos are emit- 
ted from a place of relatively high density, where they are 
matter suppressed and essentially coincident with matter 
eigenstates. As the neutrinos proceed to the outer edges 
of the star, which is at relatively low density, they pass 
through two potential resonances regions and emerge as 
mass eigenstates, see e.g. 0. The lower density res- 
onance, at ~ 10gcm~ 3 , is governed by <5mi 2 , and #12 
both of which have been measured by solar and reactor 
neutrino experiments. The higher density transition, at 
~ 10 3 g cm -3 is governed by 617113 and #13 . The for- 
mer parameter is known except for a sign (leaving the 
neutrino mass hierarchy unknown) while the later is con- 
strained by reactor neutrino experiments sin 2 613 < 0.1 
|19|. There is some uncertainty regarding how the neu- 
trinos will behave in these regions due to the unknown 
hierarchy, #13 and the density profile of the star. 

The results of various oscillation scenarios are shown in 
Table [I] where we show the ratio of the total energy flux 
in all ^s to that in v e , v e and the other four neutrinos as 
seen at a detector. A total energy flux could be inferred 
from a signal in the neutral current channel, while the 
v e and D e fluxes could be inferred through measurements 
in the charged current channels. It can be seen from the 
table that if neutrinos come from a PNS, then the energy 
flux ratios are relatively similar. The ratios for the AD- 
BH cases, on the other hand, vary considerably. If v e s 
or D e s encounter the higher resonance adiabatically, then 
there is a clear signature for the AD-BH scenario because 
the flux ratios differ by an order of magnitude from the 
PNS. In almost all other cases, there are still differences 
between the PNS and AD-BH flux ratios of a factor of 
two. 

To determine the relative fluxes in the different chan- 
nels, sufficient counts are required. In Fig. |3| we com- 
pare events that would occur in SuperKamiokande, for 
a supernova at 10 kpc, which collapses either to a PNS 
or to an AD-BH. Estimates of total counts for Super- 
Kamiokande, KamLAND, SNO (as a model for a v e sen- 
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FIG. 2: Shows the neutrino spectra for the PNS and AD-BH, in the absence of oscillations, at a distance of 10 kpc from the 
center of the object. The AD-BH neutrino spectra were calculated by summing over the entire surface of the disk including 
both trapped and untrapped regions, as described in [lH Il6|. The spectral parameters for the PNS were taken from [r?|: 
T Vc = 2.6, rj = 3.0, T„ e = 4.0, r\ = 2.8, and T„ „ T = 5.0, rj = 1.8, with an energy partition of L Ve : L Sc : L„ e of 1 : 1.3 : 7.7. 



sitive detector) and OMNIS are given in Table ITT| The 
prospects for making the distinction will depend on which 
detectors are on line at the time. It is essential to have a 
good neutral current signal measurement, as this is how 
the total flux will be measured. It is also vastly prefer- 
able to have a way to measure both the v e and v e flux if 
one wishes to distinguish between a PNS and an AD-BH. 

We have discussed the observational consequences of 
the next Galactic core collapse event forming at its cen- 
ter an accretion disk surrounding a black hole instead of 
a proto-neutron star. This type of event will produce 



significant event rates in currently on-line neutrino de- 
tectors, with a time profile that will likely differ from 
that of the canonical PNS. With a neutral current mea- 
surement, it may be possible to determine the origin of 
core-collapse supernova neutrinos by inferring the ratio 
of the total energy flux to charged current flux. For the 
AD-BH scenario this could be as large 1600:1. 
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FIG. 3: Comparison between the counts for P e + p — * e + + n at SuperKamiokande as a function of positron energy for the two 
oscillation scenarios. Plotted are the signals from the proto-neutron star and accretion disks of 0.1 M©/ s to 1.0 M©/s. With 
a normal hierarchy (left panel), a fraction of the neutrinos originally emitted as v e s are measured as £> M and v T s at the detector. 
With an inverted hierarchy and an adiabatic Q\z resonance (right panel), almost all the neutrinos originally emitted as i> e s are 
measured as i/ M and z> T s at the detector. For the PNS this makes relatively little difference as originally emitted Vp and v T s are 
also converted to v e . 
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TABLE I: Ratios of energy fluxes for various accretion disk models and oscillation scenarios as compared with the PNS. The 
PNS is calculated for the case of equipartition of energy between neutrino species, and for the luminosity ratios ratios of Fig. 
[5] The quantity v x refers to the sum of the z/ M , v T and and v T luminosities. Five oscillation scenarios are shown. H means 
that either neutrinos or antineutrinos, depending on the hierarchy, proceed adiabatically through the higher density resonance 
(i.e. Phop = 0). LA means that neutrinos proceed adiabatically through the lower density resonance, while LNA means that 
neutrino pass nonadiabatically through the lower density resonance (i.e. Phop = 1). 
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TABLE II: Integrated counts in various detectors. The columns represent the same oscillation scenarios as in Table All 
accretion disk rates may be scaled by (M/Mq) where M is the mass of material processed by the disk. OMNIS rates are quoted 
for lkt of lead perchlorate only. More detailed discussion of count rates for the PNS are in, e.g. [Til I20II21I |22| . 



